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Abstract 

Two non-directly interacting qubits with equal frequencies can become entangled via 
a Markovian, dissipative dynamics through the action of a weakly coupled Ohmic heat 
bath. In the standard weak-coupling limit derivation, this purely dissipative effect 
disappears if the frequencies are different because of the "ergodic average" used by this 
approach. However, there are physical situations where this technique is too rough to 
capture all the relevant aspects of the dissipative dynamics. In these cases, in order to 
better describe the physical behavior of the open system, it is necessary to go beyond 
the "ergodic average" . We show that, in this more refined framework, the entanglement 
capability of the environment persists also in the case of different frequencies. 

Quantum systems are usually treated as isolated: this is just an approximation, justified 
for vanishingly small couplings with the external environment. When the interaction with the 
environment is weak but not negligible, a reduced dynamics can be obtained by eliminating 
the environment degrees of freedom and by subsequently performing a so-called Markovian 
approximation pQ- [SJ. These systems are known as open quantum systems and their reduced 
dynamics is irreversible and satisfies a forward-in-time composition law: it is described 
by a so-called quantum dynamical semigroup that incorporates the dissipative and noisy 
effects due to the environment. The latter acts as a source of decoherence: in general, 
the corresponding reduced dynamics irreversibly transforms pure states (one-dimensional 
projections) into statistical mixtures (density matrices). 

One of the most intriguing aspects of quantum coherence is entanglement [6], that is the 
existence of purely quantum mechanical correlations, which has become a central topic in 
quantum information for its many applications as a physical resource enabling otherwise 
impossible information processing protocols. With reference to the entanglement content 
of a state of two qubits embedded in a same heat bath, it is generally expected that it 
would be depleted by decoherence effects. However, this is not the only possibility: if suit- 
ably engineered, the environment can entangle an initial separable state of two dynamically 
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independent systems; the reason is that, although not directly interacting between them- 
selves, there can be an environment mediated generation of quantum correlations between 
two systems immersed in it. 

This possibility has been demonstrated analytically for two qubits with a same oscillation 
frequency [8], d [9] and two identical harmonic oscillators [10] evolving according to a reduced 
master equation of the typical Lindblad form [TU [T2] . obtained via the so-called weak- 
coupling limit. This technique is based on the fact that the dissipative effects are visible 
only over a coarse-grained time-scale At, so large that the free dynamics of the embedded 
systems can be averaged out over At [2]- 03- The resulting elimination of too rapid oscillations 
is mathematically implemented through a time "ergodic average" . 

Noticeably, such a prescription guarantees that, unlike for reduced dynamics of Redfield 
type (see [HIE]), the resulting quantum dynamical semigroups consist of completely positive 
maps [H [21 El ED] ■ Complete positivity ensures that the open quantum evolution is consistent 
with entanglement in the sense that not only the positivity of any initial density matrix of the 
open system is preserved in time, but also that of any initial state of the open system coupled 
to any other possible ancillary system. Indeed, only complete positivity can guarantee the 
full physical consistency of the Markovian approximations that one uses to describe an open 
quantum dynamics; in other words, without complete positivity, it always occurs that at 
least one initial state carrying entanglement between the given system and an ancilla will 
assume negative eigenvalues in the course of time [5]. 

It turns out that, in the standard weak-coupling limit approximation, when two qubits 
or two harmonic oscillators embedded in a same environment have different oscillation fre- 
quencies lo\ 7^ u>2, no matter how small the difference u>2 — oji is, the elimination of too rapid 
oscillations destroys the generation capability that the environment possesses when oj\ = UJ2- 

In this paper, we study this behavior in the case of two qubits weakly interacting, via a 
Ohmic coupling, with a heat bath made of free bosons at high temperature. We shall first 
relate the sharp dependence of the entanglement capability of the environment to the oscil- 
lation frequencies of the two qubits on the drastic elimination of too fast oscillations through 
the "ergodic average" . This procedure is only allowed when the coupling to the environment 
is such that the coarse-grained time-scale At can effectively be considered infinitely large. 
However, there are situations where this approximation is not physically meaningful and At 
should be kept finite. In these cases, to better capture the effects of the open system dynam- 
ics, a less rough time coarse-graining is needed; in the following, we use a master equation 
derived without recourse to the ergodic average that nevertheless generates a completely 
positive dynamical semigroup. We show that, in this refined framework, the entanglement 
generation capability of the environment is preserved even when u\ 7^ u>2- 

The problem we will address in the following regards whether two non-interacting qubits 

1 For the system composed by two harmonic oscillators immersed in a heat bath of other oscillators similar 
results have been obtained by studying numerically their exact time evolution [T3] . 
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with Hamiltonian 



Wl (1) 
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can become entangled when weakly coupled to free spinless Bosons in thermal equilibrium 
via a (finite volume) interaction of the form 



5/ = AX(/)(af ) + a®) , X(f) = J2(f(k)al + r(k)a k 

k 



(2) 



where a k and a k denote the creation and annihilation operators of Bose modes with momen- 
tum k and energy u(k), f(k) is a one-particle Bose state in momentum representation, A is 
a small, dimensionless coupling constant, while a^l = ® 1 and a^l — 1 ® 01,3 represent 
the first and third Pauli matrices for the two qubits. d The total system comprising the two 
qubits and the thermal bosons will thus be described by the Hamiltonian 



H = H s + H B + Hi , H B = ^( k ) 4 a k- 



(3) 



In general, the state of the compound system S + B at time t > is a correlated state 
Pse(t) from which the state of the two qubits can be extracted as p(t) = Tr B(psB(t)) by 
tracing out the environment degrees of freedom; however, the time-evolution equation for 
p(t) is quite complicated. In order to arrive at a memoryless master equation, one starts with 
the physically acceptable hypothesis that the initial state of open system plus environment 
be of the uncorrelated form p®p/3, where pp oc exp {—(3Hb) is the Boson thermal equilibrium 
state at inverse temperature f3. Then, one goes to the interaction representation by replacing 
pit) with psB{t) = e lt ( Hs+HB ^ psB{t) e~ lt ( H s+ H B) anc i l OQ ks at the time-evolution over time- 
intervals At by stopping at the first significative order in A in the time-ordered expansion of 
p{t + At). Then, 



pit + At) - pit) 
At 



A_ 2 
At 



t+At 



dtoTr 



Hjih), Hi(t 2 ),p SB {t) 



where 

H I (t)=X t (f) (aP(t)+a?\t) 



(4) 
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with cr[ (t) = cosiuj a t) a[ a> — sin(cj a t) a. 

Notice that the variation of p(t) is of order A 2 , namely it is non-negligible only over times 
t = t/A 2 ; if A 1 and At is such that A 2 At is small on the scale of r, one may reasonably 
substitute in (HI) the finite ratio with a time-derivative: 



(a) 
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1,2. 



p((r + A 2 At)/A 2 )-p(r/A 2 ) 



At 



dtPit) 



(6) 



2 For sake of simplicity, we consider the two qubits located at a same point in space. For the effects on 
entanglement creation that result when the two qubits are spatially separated we refer to E] for equal 
frequencies and to p~5] in the case of u)\ ^ u>2 ■ 
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Indeed, the error is bounded by A 2 At. Further, if At is chosen much larger than the decay 
time tb of the environment two-point time-correlation functions, one may approximate psB{t) 
with the uncorrelated state pit) £g> p@ as it was at t — [21 Ej- Therefore, if tb <C At, one 
gets the following approximated master equation (in interaction representation): 



^2 ft+At / \ 

dtpit) = J dtx J dt 2 Trf [^/(ti), [Hj(t 2 ), pit) ® Pp\\ J • (7) 

The meaning of At is that of a time-coarse graining parameter naturally associated to 
the slow dissipative time-scale r = t\ 2 . More precisely, significant variations of the system 
density matrix due to the presence of the environment can only be seen after a time At = 
A 2 At has elapsed. Given the coupling constant A 1, the dissipative time-scale is set and 
actual experiments cannot access faster time-scales; furthermore, if At C 1 on the scale 
of r, then the experimental evidences are consistently described by a master equation as 
in (J7J). The weak-coupling limit consists in letting A — > [21 S]; then, variations of the 
system density matrix as in ([7]) are actually visible only if At — > +oo. This allows one to 
eliminate all fast oscillations through a time "ergodic average" as in the standard weak- 
coupling limit approach [HIE]- Instead, if the system-environment coupling A is small, but 
not vanishingly small, then At cannot be taken infinitely large. In such cases, the usual 
weak-coupling limit techniques provide an approximation which is too rough to properly 
describe the dissipative time-evolution and one needs a more refined approach in order to be 
able to keep contributions that would otherwise be washed out. El 

Returning to the Schrodinger representation, (J7|) yields the following memoryless master 
equation (see [TU (TTJ [18] ) of Kossakowski-Lindblad type [21 E]: 

d t p(t) = -i [h s + A 2 H At , p(tj\ + D At [p(t)] . (8) 

The environment contributes to the generator of the reduced dynamics with an Hamiltonian 
H At and a purely dissipative term D At [p(t)}] both of them depend on the environment 
through the two-point time-correlation functions 

G(t) = Tr( P/3 X(/)X,(/)) = G(-t)* . (9) 

The bath-induced Hamiltonian H At contains a bath-mediated qubit-qubit interaction 

H% = £ h tJ (At)a^>af , (10) 

i,j=l,2 



3 A different approach is developed in |16j : there, a non-Markovian weak coupling approximation of 
the reduced dynamics is introduced, leading to a two-parameter family of dynamical maps, with a time- 
dependent generator [2] . We stress that this treatment is completely different from the one discussed below, 
which instead describes the reduced two-atom dynamics in terms of a Markovian, one parameter semigroup. 
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where the 2x2 matrix h&t — [hij(At)] is real, but not necessarily Hermitian. It is convenient 
to introduce the following matrices [& e j] = 2 ^i ^ ij ' £ = then 

h At = *t H { £ ] V + (V #gi) ^ T } (11) 
where = [H^, (At)], a, b = 1,2, is explicitly given by 

fAt pAt 

H^(At) = -— J dt a j( dt 2 e^'-^sign^ - - h) , (12) 

where X T denotes matrix transposition. Instead, the purely dissipative part can be written 



as 



o,6=l,2; i,i'=l,2 V 



4 ] ^\p(t)\\ , (13) 



where the 2x2 matrices = [C^(At)] read 

Cff = ^ff* , (14) 



in terms of the matrices -D^t = [-^w (AO] with entries 

i pAt pAt 

D^(At) = — dt x dhe^^'^Gih-h) . (15) 

At Jq Jq 

The 4x4 Kossakowski matrix C&t formed by the 2x2 blocks a,b = 1,2, involves the 

Fourier transforms of the correlation functions ([9]) and is automatically positive definite; this 
fact guarantees that the master equation (jBJ) generates a semigroup of dynamical maps 7 t A * 
on the two-quibit density matrices which are completely positive PQ-[S], whence, as discussed 
in the introduction, fully physically consistent. 

Given two qubits weakly interacting with their environment, a sufficient condition for 
them to get entangled at small times by the completely positive reduced dynamics has been 
derived in [TUJ [20] and is based on the properties of the generator in (jSJ), that is on the 
interaction Hamiltonian (1101) and the dissipative contribution (1131) . We shall focus on an 
initial two qubit state of the form | J.) ® | |), where | I), | |) are the eigenstates of cr 3 ; then, 
the condition is as follows: 
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5 = D {u J(At) D { +1 ] (At) - V ( £> + iH^ < , (16) 
where, from (|T2|) and (|T5l) . 

(12) D^(At) + (D { ^(At)Y 
u At — 2 ^ ' 

= H^iAt) +H ( +l\At) . (18) 



,(12) .-,(12) 
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If 5 > 0, such a separable state cannot get entangled by the environment at small times [20]. 

We are interested in the capacity of the environment to generate entanglement, at least 
at small times with respect to the dissipative time-scale. Because of the positivity of the 

= 1,2, are positive 



Kossakowski matrix C^t the diagonal block matrices C^ t and D 



'At 



definite; therefore, the dissipative entanglement generation depends on the quantity T> 
but also on how it interferes with in ffTUj) . In the case at hand, it turns out that 



(12) 
At J 



W12) -1/(12) 



l At 



V 



(12) 
At 



+ 



H 



(12) 
At 



(19) 



Indeed, using ((9]), one can check by explicit computation that the quantities e Wl ) A */ 2 j) 



and e Ul ^ A ^ 2 Ti.^' are both real. We shall thus concentrate on the purely dissipative 



(12) 
At 



entanglement generation, that is on the difference 



D_ 



V 



(12) 
At 



(20) 



If 6 is negative, a fortiori also 5 in (1161) is negative. 

We shall consider (infinite volume) Ohmic correlation functions 



G(t) = J du) e~ u) / ulc io ^coth --- cos tut — i sin cut J , 



(21) 



where u c is a Debye cut-off; then, setting sinc(x) 
explicitly gets 



SIM 

X 



and (a, e) = (1, — ), (2, +), one 



eAt 



-\u>\/ul c 



smc 



(uj - Lo a )At 



V 



(12) 
At 



■y / dujue~ H/uJc coth(/9u;/2)ffliK- 



(uj — (jJi)At 



smc 



(uj — uj 2 )At 



(22) 



(23) 



The weak-coupling limit amounts to At — ► +oo; since lim —smc(a(x — x )) = 5(x — x ), 
one finds 



lim D^(At) = 4tt£^- 



UJa/uic 



lim 

At— >+oo 



V 



(12) 
At 



27r5 WlWa uj 2 e-" 2/Wc coth 



/3oj 2 



(24) 
(25) 



If u>i = uj 2 , the difference ( J201) is always negative: 5 = — An 2 uj 2 e 2w2 / Wc ; instead, if c^i 7^ cj 2 , 
5 = (At) (At) > 0. Thus, in the latter case, the initial separable state | |) <8> | T) 
cannot get dissipatively entangled by the environment at small times, unless u>\ = uj 2 . 
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This behavior is characteristic of physical situations where an infinitely large coarse- 
grained time-scale At is justified by a vanishingly small coupling A to the environment. 
However, if A is small, but not negligibly small, the qubit density matrix effectively varies 
over large but finite At. Then, terms of order 1/At like 



T)(12) 



sinc(<5u;At) u a e -" a/ " c coth(/3u; a /2) (26) 



7T 

a=l 



can make 5 negative even when 5u = (u 2 — ^i)/2 > 0. Indeed, for high temperatures 
/3&i,2 "C 1 and large cut-offs toi^/^c "C 1, expanding (12^1) and (1261) yields @ 



16vr 2 



2 



(l - (35lo - sinc 2 (^At)) . (27) 



If the qubit frequency difference 5u and the coupling strength A are such that 5uAt <C 1 
then a further expansion of the sine function yields 

i6. V ^_(MW ); (28) 



P 

, . i . . 1 n su(Aty 

which is negative when ^> p > 



uj 1j2 3 

In conclusion, the fact that the environment may generate two-qubit entanglement only if 
the qubit frequencies are equal, is a consequence of the weak-coupling limit and the associated 
time "ergodic average"; in such a case, the coupling constant A — > so that the coarse-grained 
time-scale At over which the qubit density matrix effectively changes due to the presence of 
the environment become so large that all off- resonant phenomena are averaged out. Instead, 
if the physical conditions ask for a coupling constant A that is small but not vanishingly so, 
a consistent description of the open dynamics requires At finite; in this way, one may keep 
track of finer effects and save the possibility of dissipative entanglement generation even if 
0J\ 7^ u 2 . 



4 The high temperature hypothesis has been made in order to permit an analytical study of the behavior 
of 5. This is the worst possible scenario; indeed, in general, lower temperatures favor the entanglement 
creation capability of the environment [15] . 
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